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ABSTRACT 



Context. The general properties (e.g., activity class, star formation rates, metallicities, extinctions, average ages, etc) of luminous and 
ultraluminous infrared galaxies (LIRGs and ULIRGs, respectively) in the local universe are well known since large samples of these 
objects have been the subject of numerous spectroscopic works over the last three decades. There are, however, relatively few studies 
of the spatially-resolved spectroscopic properties of large samples of LIRGs and ULIRGs using integral field spectroscopy (IFS). 
Aims. We are carrying out an IFS survey of local (z < 0.26) samples of LIRGs and ULIRGs to characterize their two-dimensional 
spectroscopic properties. The main goal of this paper is to study the spatially resolved properties of the stellar populations and the 
excitation conditions in a sample of LIRGs. 

Methods. We analyze optical (3800-7200A) IFS data taken with the Potsdam Multi-Aperture Spectrophotometer (PMAS) of the 
central few kiloparsecs of eleven LIRGs. To study the stellar populations we fit the optical stellar continuum and the hydrogen 
recombination lines of selected regions in the galaxies. We analyze the excitation conditions of the gas using the spatially resolved 
properties of the brightest optical emission lines. We complemented the PMAS observations with existing //ST/NICMOS near-infrared 
, continuum and Pact imaging. 

Results. The optical continua of selected regions in our LIRGs are well fitted with a combination of an evolved (~ 0.7 - 10 Gyr) stellar 
population with an ionizing stellar population (1 - 20Myr). The latter population is more obscured than the evolved population, and 
has visual extinctions in good agreement with those obtained from the B aimer decrement. Except for NGC 7771, we found no clear 
evidence for an important contribution to the optical light from an intermediate-aged stellar population (~ 100 - 500 Myr). Even after 
correcting for the presence of stellar absorption, a large fraction of spaxels with low observed equivalent widths of Ha in emission 
still show enhanced [Nii]/16584/Hq' and [S ii]/1/167 17,673 1/Hor ratios. These ratios are likely to be produced by a combination of 
photoionization in H n regions and diffuse emission. These regions of enhanced line ratios are generally coincident with low surface 
brightness H n regions and diffuse emission detected in the Ho- and Paa images. Using the PMAS spatially resolved line ratios and 
the NICMOS Pao- photometry of H n regions we find that the fraction of diffuse emission in LIRGs varies from galaxy to galaxy, and 
it is generally less than 60% as found in other starburst galaxies. 
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1. Introduction 2002; Sanders & Ishida 2004; Scoville et al. 2000; Bushouse 
«J ■ et al. 2002; Alonso-Herrero et al. 2006), the nuclear activ- 
<3 . Luminous and Ultraluminous Infrared Galaxies (LIRGs and hy dass> star formation rates> and extinctions (e.g., Heckman, 
ULIRGs^ respectively) with infrared 8 - 1000 /an luminosities AmuS; & MOey 1987; Amus> H eckman, & Miley 1989; Kim 
L IR -10 -10 L G andL IR -10 -10 L Q , respectively (see et al 1998; Veilleux et al 1995i 1999; Qoldader et 
Sanders & Mirabel 1996) are among the most luminous objects al 1995; Wu et al 1998; Rigopoulou et al . 1999; Heckman 
in the local universe, and are believed to be powered by strong et al 2000; Yuan> Kewley ^ & Sanderfj 2010X steUar popula . 
star formation and/or AGN activity. The main properties of local tionfj ( p oggiant i & W u 2000; Chen et al. 2009; Rodriguez 
LIRGs and ULIRGs are well known since these two classes of Zmdn ^ Tadhuntel - & Gonza i ez Delgado 2009, 2010a), metal- 
galaxies have been extensively studied using imaging and spec- licitiefj ( Rupk ^ Veilleux & Baker 20 08), and mo i eC ular gas 
troscopy over the last three decades. These properties include, CQntent ( Sanders et al 1991; Gao & Solomon 2004; Gracia- 
among others, morphologies (e.g., Veilleux, Kim, & Sanders Carpicj et al 2006 2008) 

* Based on observations collected at the German-Spanish There are. however, relatively few studies of the spatially- 
Astronomical Center, Calar Alto, jointly operated by the Max- resolved spectroscopic properties of LIRGs and ULIRGs us- 
Planck-Institut fur Astronomie Heidelberg and the Instituto de ing integral field spectroscopy (IFS). Most works have focused 
Astrofisica de Andalucfa (CSIC). on individual famous objects or small samples (e.g., Colina, 
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Table 1. Statistical properties of the spatially resolved of line emission 
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0.65 


0.74 


0.44 
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0.38 


207 


0.09 


0.15 


0.15 
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20. 


214 
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81 
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UGC 1845 


1.16 


1.24 
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20 
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1.39 


0.84 


171 


0.39 
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0.37 
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0.07 
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11. 


11. 


178 


NGC 2388 


0.52 


0.64 


0.61 


79 


0.58 


0.68 


0.37 


171 


0.42 


0.51 


0.44 


150 


0.06 


0.13 


0.16 


79 


16. 


22. 


17. 


171 


MCG +02-20-003 


0.66 


0.82 


0.47 
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0.44 


0.51 


0.37 
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0.44 
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0.12 
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22. 
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23. 


220 


NGC 5936 


0.39 


0.45 


0.26 


71 
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0.44 


0.09 
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0.30 


0.32 


0.08 
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0.05 
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16. 


256 


NGC 6701 


0.58 


0.57 


0.20 


28 
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0.80 


0.46 


250 


0.45 


0.61 


0.57 


177 


0.08 


0.15 


0.19 


46 


9. 


12. 


10. 


250 


NGC 7469* 


1.17 


1.76 


1.59 


77 


0.64 


0.72 


0.37 


235 


0.51 


0.57 


0.24 


196 


0.07 


0.09 


0.09 


75 


20. 


27. 


19. 


234 


NGC 7591 


0.58 


0.92 


0.98 


29 


0.54 


0.70 


0.45 


241 


1.03 


1.29 


0.86 


208 


0.19 


0.26 


0.21 


70 


7. 


11. 


11. 


242 


NGC 7771 


0.59 


0.73 


0.50 


81 


0.70 


0.90 


0.65 


270 


0.74 


0.89 


0.67 


254 


0.13 


0.19 


0.19 


125 


7. 


13. 


17. 


272 


All 


0.66 


0.85 


0.77 


739 


0.56 


0.71 


0.51 


2261 


0.52 


0.71 


0.64 


1852 


0.08 


0.14 


0.15 


811 


17 


24 


24 


2261 



Notes. — For each line ratio and the EWtHa)^ we list the values of the median, average, standard deviation and number of spaxels with measurements. The statistics is done using the 
values of the individual spaxels and thus are not light weighted. That is, the average value of a given line ratio is not necessarily equal to that measured from the integrated spectra of the 
galaxy. The line ratios are not corrected for extinction or stellar absorption. 
The statistics is only for the narrow component of the hydrogen lines. 



Arribas, & Borne 1999; Arribas, Colina, & Clements 2001; 
Murphy et al. 2001; Li'pari et al. 2004a,b; Colina, Arribas, & 
Monreal-Ibero 2005; Monreal-Ibero, Arribas, & Colina 2006; 
Garcfa-Marin et al. 2006; Reunanen, Tacconi-Garman, & Ivanov 
2007; Bedregal et al. 2009). IFS instruments working in the op- 
tical and infrared spectral ranges on 4 and 8 m-class telescopes 
are now ubiquitous, and thus this situation is changing rapidly. 
Our group is involved in an ambitious project intended to charac- 
terize the detailed optical and infrared spectroscopic properties 
of local samples of LIRGs and ULIRGs using a variety of IFS 
facilities. 

This is the second paper in a series presenting obser- 
vations with the Potsdam Multi-Aperture Spectrophotometer 
(PMAS, Roth et al. 2005) of the northern portion of a sam- 
ple of local LIRGs defined by Alonso-Herrero et al. (2006). 
The PMAS sample is in turn part of the larger IFS survey 
of nearby (z < 0.26) LIRGs and ULIRGs in the northern 
and the southern hemispheres using different IFS instruments. 
The optical ones include, apart from the PMAS instrument, the 
VIMOS instrument (Le Fevre et al. 2003) on the VLT, and the 
INTEGRAL+WYFFOS system (Arribas et al. 1998; Bingham 
et al. 1994) on the William Herschel Telescope (WHT). The 
first results of this project can be found in Alonso-Herrero et 
al. (2009, hereafter Paper I) for the PMAS atlas of LIRGs, 
Garcia-Marfn et al. (2009a,b) for the INTEGRAL results of 
ULIRGs, and Arribas et al. (2008), Monreal-Ibero et al. (2010a) 
and Rodriguez Zaurin et al. (2010b) for the VIMOS results. 
Additionally, Pereira-Santaella et al. (2010) studied the spatially 
resolved mid-infrared properties of LIRGs using the spectral 
mapping capability of IRS on Spitzer, and Bedregal et al. (2009) 
presented a detailed near-IR IFS study of a local LIRG using 
SINFONI on the VLT. 

In this paper we study in detail the stellar populations, ex- 
citation conditions, and diffuse emission in the central regions 
(a few kiloparsecs) of a sample of 1 1 LIRGs using the PMAS 
data. The paper is organized as follows. We present the obser- 
vations and data analysis in Section 2. In Section 3 we describe 
the modelling of the stellar populations. The results regarding 
the morphology, stellar populations, excitation conditions, and 
diffuse emission are presented in Sections 4, 5, and 6. We give 
our conclusions in Section 7. 

2. Observations and Data Analysis 

2.1. Sample, observations and Data Reduction 

The observations and data reduction of the PMAS data are de- 
scribed in detail in Paper I. Briefly, we used PMAS on the 



3.5 m telescope at the German-Spanish Observatory of Calar 
Alto (Spain) to observe a sample of 11 local LIRGs. These 
comprise the majority of the northern hemisphere portion of the 
volume-limited (v = 2750 - 5200kms _1 ) sample of LIRGs of 
Alonso-Herrero et al. (2006). This volume-limited sample was 
originally drawn from the IRAS Revised Bright Galaxy Sample 
(RBGS, Sanders et al. 2003). The range of IR luminosities of the 
PMAS sample is log(L IR /L ) = 11.05 - 11.59, and the galaxies 
are at an average distance of 61 Mpc (for the assumed cosmol- 
ogy Ho = 75kms Mpc -1 ). As the sample is flux-limited, it 
is composed mostly by moderate IR luminosity systems, with 
an average log(Li R /L G ) = 11.32 for the full sample of Alonso- 
Herrero et al. (2006). 

The PMAS observations were taken with the Lens Array 
Mode configuration, which is made of a 16 x 16 array of mi- 
crolenses coupled with fibers called hereafter spaxels. We used 
the 1" magnification that provides a field of view (FoV) of 
16" x 16". All the galaxies were observed with a single point- 
ing, except NGC 7771 for which we obtained two pointings to 
construct a mosaic with a 28" x 16" FoV. We covered the wave- 
length range of 3800 - 7200 A, using the V300 grating with a 
spectral resolution of 6.8 A full width half maximum (FWHM). 
The full description of the data reduction procedure can be found 
in Paper I. 

In addition to the PMAS data, in this paper we make use 
of the HST/NICMOS Fl 10W (Ac = 1.1 /im) and F160W (A c = 

1 .6 /im) continuum observations to construct continuum near-IR 
color maps, as discussed by Alonso-Herrero et al. (2006). The 
only additional step needed for the reduced NICMOS images 
was to rotate and trim them to match the orientation and FoV, 
respectively, of the PMAS images. 



Table 2. Definition of indices of the stellar absorption features. 



Name 


Blue Cont. 
Ac AA 


Red Cont. 

A, AA 


Line 

/t c AA 


Ref 


C6500/C4800 


6500 


50 


4800 


50 






1 


D„(4000) 


3900 


100 


4050 


100 






2 


H<5 A 


4060 


38.5 


4145 


32.5 


4102 


38.5 


3 


Wa 


4776 


12 


4948 


12 


4861 


30 


4 


Ha A 


6510 


8 


6616 


8 


6563 


30 


4 



Notes. — All the wavelengths are in A. 

References: 1. Kim et al. (1995) and Veilleux et al. (1995). 2. Balogh et 
al. (1999). 3. Worthey & Ottaviani (1997). 4. Gonzalez Delgado et al. 
(2005). 
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2.2. PMAS spatially resolved emission line ratios and 
equivalent widths 

We constructed spectral maps of the brightest emission lines 
in an automated fashion using our own IDL routines as 
well as the IDL-based MPFITEXPR algorithm^] the lines to 
Gaussian functions and the adjacent continuum to a straight 
line, on a spaxel-by-spaxel basis, as described in more detail 
Paper I. These fits provided the flux, equivalent width (EW), 
and full width half maximum of the emission lines. In this 
paper, we generated spectral maps of the following optical 
line ratios: [Ora]^5007/H/3, [Oi]i6300/Ha, [Nn],l6584/Ha, and 
[Sn]/U6717,6731/Ha. We also generated maps of the EW of 
the Ha line in emission (EW(Ha) em ), which are thus measured 
as positive numbers. The maps of the line ratios and EW are 
not corrected for the presence of Ha and H/3 in absorption. The 
spectral maps of the observed line ratios together with those of 
the observed flux of Ha and EW(Ha) em are shown in Figure 1, 
except for IC 860, for which the line emission is compact (see 
Paper I). In Table 1 we list for each of the emission line ratios 
the number of spaxels where it was possible to obtain a measure- 
ment, the median, the average value and the standard deviation, 
for each galaxy and the full PMAS sample. 

We also produced standard optical diagnostic diagrams us- 
ing the brightest emission lines (Baldwin, Phillips, & Terlevich 
1981; Veilleux & Osterbrock 1987) on a spaxel-by-spaxel ba- 
sis for the galaxies in our sample. These diagrams provide use- 
ful information on the excitation conditions of different regions 
in galaxies, such as, photoionization by young stars, shocks, 
and AGN photoionization. In Paper I we presented such dia- 
grams for the nuclear and integrated measurements of our sam- 
ple. The spatially-resolved diagnostic diagrams for each of the 
galaxies in our sample are shown in Figure 2. We also produced 
[Nn]^6584/Ha vs. [S n],U6717,6731/Ha diagrams (Figure 3) 
for each of the galaxies in our sample. These diagrams have an 
advantage over the Veilleux & Osterbrock (1987) diagrams in 
that they contain more data points, as the [O rn]/l5007/Hj6 line ra- 
tio can be strongly affected by both the presence of an underlying 
stellar absorption and extinction. As we shall see in Sections 5 
and 6, the effects of the Balmer absorption stellar features on the 
observed line ratios are not negligible in regions with low val- 
ues of EW(Ha) em or EW(Hy6) em . For this reason, the individual 
measurements in diagrams of Figures 2 and 3 are color-coded 
according to arbitrarily chosen ranges of EW(Ha) em . All these 
diagrams and the effects of the correction for stellar absorption 
features will be discussed in Section 6. 

The errors of the line ratios depend on the observed val- 
ues of EW(Ha) em and the S/N of the spectra. We estimated the 
typical uncertainties by comparing the line ratios measured au- 
tomatically with our IDL routines with those fitted manually 
with the splot routine within IRAF for selected spaxels in each 
galaxy. For each galaxy, the comparison was made for spaxels 
within the smallest observed range of EW(Ha) em where the un- 
certainties are the highest. By choosing spaxels with low val- 
ues of EW(Ha) em , we basically estimated an upper limit to the 
uncertainties of the observed line ratios. As can be seen from 
Figure 3 the largest uncertainties in the [Nii]/16584/Ha and 
[Sn]X16717,6731/Ha ratios are 10 - 25%, and 15 - 40%, re- 
spectively, depending on the range of EW(Hff) em and the galaxy. 



1 http://www.purl.com/net/mpfit developed by Markwardt (2008). 
We fitted 



2.3. Extraction and Analysis of the 1D spectra of selected 
regions 

For each galaxy we extracted the nuclear and integrated spec- 
tra as done in Paper I. Briefly, we identified the position of 
the optical nucleus as the peak of the 6200 A continuum emis- 
sion, and extracted the nuclear spectrum using the corresponding 
spaxel. The physical size covered by the nuclear spectrum for 
each galaxy was given in Paper I, and it is typically the approxi- 
mate central 300 pc. The integrated spectrum of each galaxy was 
extracted by defining ~ 6200 A continuum isophotes and then 
summing up all the spaxels contained within the chosen external 
continuum isophote to include as much as possible of the PMAS 
FoV. The area covered by the integrated spectrum can be seen 
in figure 1 in Paper I, and it is generally the central 5-8 kpc, 
depending on the galaxy. 

For the three galaxies observed under photometric condi- 
tions: NGC 23, NGC 2388, and NGC 7771, we extracted spectra 
of a number of bright Hn regions. The locations of the selected 
regions are shown in the lower panels of Figure 4. These regions 
were chosen to probe, for each galaxy, a range of EW(Ha) em 
and continuum slopes as can be seen from Figure 5. Veilleux 
et al. (1995) found that the extinction in (U)LIRGs powered by 
star formation is correlated with the shape of the optical contin- 
uum. We estimated the shape of the optical continuum by mea- 
suring the ratio (C6500/C4800) between the continuum fluxes (see 
Figure 4) near H/3 and Ha at 4800 and 6500 A, respectively. The 
definition of this ratio is given in Table 2. Similarly the gas ex- 
tinctions are related to the HST/NICMOS F160W - Fl 10W col- 
ors (see Alonso-Herrero et al. 2006). Figure 4 shows how the 
selected Hn regions in these galaxies span a range in F160W to 
Fl 10W flux ratios, and there is a good correspondence between 
red near-IR colors and steep optical continua. Table 3 gives de- 
tails of the observed properties of these regions. Given the good 
quality of the spectra for these three galaxies, we will use them to 
do a detailed modelling of the stellar populations in these galax- 
ies (Section 3.1). 

As we will explain in Section 5.1, we can use the 4000A- 
break and the Balmer lines in absorption, and in particular the 
H6 line, to infer an estimate of the average age of the stellar 
populations. Although the H6 nebular emission line is much 
weaker than Ho- and H/3, it can also be observed in emission 
in bright/young H 11 regions (see for instance the HII-1 region in 
NGC 23, Figures 4 and 5). To minimize the contamination from 
the Balmer lines in emission, for each galaxy we selected regions 
with low values of EW(Ha) em so we could attempt to measure 
the HS feature in absorption. We produced the average spectrum 
of this low EW emission by summing up the spaxels with the 
specified range of EW(Ha) em for each galaxy. The individual 
spaxels had typically 20 A < EW(Ha) em < 6 A, but the specific 
range depended on the galaxy and the S/N of the data. From 
the comparison between the maps of the observed EW(Ha) em in 
Figure 1 and the //5T/NICMOS Paa maps shown in Paper I, it 
is clear that the low EW(Ha) em values are associated with re- 
gions of diffuse emission or low surface brightness H n regions. 
Table 4 lists the values of EW(Ha) em as measured from the aver- 
age spectrum of regions of low EW of Ha of each galaxy. All the 
extracted spectra were shifted to rest-frame wavelengths prior to 
the analysis and fitting the stellar populations. 

In Section 5 we will study the stellar populations of LIRGs 
using the 4000 A-break and the HS stellar feature. We adopted 
the definition of Balogh et al. (1999) for the D„(4000) index and 



4 Alonso-Herrero et al.: PMAS Optical Integral Field Spectroscopy of Luminous Infrared Galaxies. 

Table 3. Observed properties of selected star-forming regions. 



Region 


Pos X 
arc sec 


Pos y 
arcsec 


Size 
pc x pc 


EW(Ha) em 
A 


D„(4000) 


H<S A 


C6500/C4800 


NGC 23 


HII-1 


-3.0 


-3.5 


289 x 578 


93.9 


1.37 




1.14 


HII-2 


-3.5 


+2.5 


578 x 578 


50.8 


1.31 


1.95 1 


1.12 


HII-3 


+0.5 


+3.5 


578 x 578 


47.2 


1.20 


2.68+ 


0.98 


HII-4 


+ 1.5 


-2.5 


578 x 578 


26.9 


1.22 


5.01 


0.80 


NGC 2388 


HII-1 


+4.5 


-0.5 


560 x 560 


66.3 


1.27 


3.87 


1.41 


HII-2 


-2.5 


+0.5 


560 x 560 


64.0 


1.35 




1.69 


NGC 7771 


HII-1 


-2.5 


-2.0 


277 x 554 


89.0 


1.18 




0.99 


HII-2 


-5.5 


-0.5 


554 x 554 


39.8 


1.46 




1.57 


HII-3 


-3.5 


+ 1.5 


554 x 554 


26.8 


1.44 




1.45 


Nucleus2 


-2.0 


+ 1.0 


277 x 277 


34.0 


1.21 


2.00 


1.48 



Notes. — The positions of the extracted regions are given relative to that of the nucleus of the galaxy. All the reported values of EW(HaO cm are in 

emission. 

The values of H5 A are not corrected for H<5 nebular emission. f The nebular H<5 line is seen in emission within the stellar absorption feature. 
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Fig. 4. Upper panels: ffiiT/NICMOS 1 .6 pm to 1 . 1 pm flux ratio maps. Dark regions show a deficit of 1 . 1 pm flux, that is, a redder 
F160W - F110W color. All three galaxies are shown on the same flux ratio scale. Lower panels: The contours are the observed 
Ha fluxes of the three galaxies for which extracted spectra of circumnuclear Hn regions. For each galaxy we show the regions of 
interest: the nucleus (defined as the peak of the 6200 A continuum emission, light blue dots), bright Hn regions (red, blue, purple, 
and green dots). The orange dots are the spaxels used to produce the average spectrum of regions with relatively low EW of Ha 
(see Section 2.3). For NGC 7771 we additionally extracted a spectrum of the possible location of the true nucleus (shown as a black 
dot), based on the Ha velocity field (see Paper I for details). 



that of Worthey & Ottaviani (1997) for the H6 A index0 We mea- 
sured these indices for the nuclear and integrated spectra, as well 
as the average spectra of regions of low EW(Ha) em . The pseudo- 
continuum bands of these two indices, as well as the line window 
for the H<5a index are listed in Table 2. For the majority of the 
selected Hn regions (see Table 3), as well as the nuclei and in- 
tegrated emission (Table 4), the measured values HSa are only 
lower limits as the observed values of EW(Ha) em in emission 



2 We note that the value of this index is positive for absorption, and 
throughout this paper we will use this symbol to indicate the feature in 
absorption. 



imply a contribution from nebular H6 emission to the index. In 
Tables 3 and 4 we marked those regions where we clearly de- 
tected the nebular H<5 line in emission within the stellar absorp- 
tion. We also measured the Ha flux and the EW(Ha) em of the 
line in emission in all the extracted spectra. 



3. Modelling of the Stellar Populations 

As discussed in Paper I the nuclear and integrated spectra of our 
sample of LIRGs show evidence for the presence of an ioniz- 
ing stellar population, plus a more evolved stellar population 
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Table 4. Observed properties for the nuclei, average of regions with low EW(Ha) em and integrated emission. 



Galaxy Nuclear Low-EW(Ha) cra Region Integrated* 

EW(Hor) em D„(4000) H<5 A EW(Ha) cra D„(4000) H<5 A EW(Hor) em D„(4000) H<5 A 



AAA 


NGC 23 


12.8 


1.40 


3.80 


7.3 


1.34 


3.54 


32.8 


1.29 


3.13+ 


MCG +12-02-001 


145.6 












92.5 






UGC 1845 


20.4 


1.60 


6.66 


10.5 






25.1 






NGC 2388 


53.1 


1.30 


2.58+ 


10.0 


1.40 


4.22 


28.8 


1.35 


3.21 


MCG +02-20-003 


143.6 


1.42 


2.01+ 








39.9 


1.36 


3.85 


IC 860 


1.5 


1.42 


5.67 










1.36 


4.40 


NGC 5936 


66.1 


1.33 


2.29 


16.0 


1.25 


2.92 


38.6 


1.12 


2.40+ 


NGC 6701 


34.8 


1.23 


2.20+ 


14.0 


1.24 


2.58 


17.2 


1.15 


2.57 


NGC 7591 


31.4 


1.55 










9.2 


1.39 


2.15 


NGC 7771 


23.0 


1.18 


5.27 


12.0 


1.28 


4.33 


15.5 


1.36 


2.22+ 



Notes. — 'Integrated spectra refers to the central ~ 5 kpc to ~ 8 kpc regions, depending on the galaxy (see Paper I for details). 
All the reported values of EW(Ha') cra are for the line observed in emission. 

The values of HSa are not corrected for H<5 nebular emission. +The nebular H<5 line is clearly seen in emission within the stellar absorption feature. 



as indicated by the presence of strong absorption stellar fea- 
tures in the blue part of the spectrum. This is also apparent for 
the H ii regions in our sample of LIRGs, including those with 
the largest EW(Ha) em (see Figure 5). These strong absorption 
stellar features appear to be a general property of local LIRGs 
and ULIRGs (see e.g., Armus et al. 1989; Veilleux et al. 1995; 
Kim et al. 1995, 1998; Marcillac et al. 2006; Chen et al. 2009; 
Rodriguez Zaurin et al. 2009), as well as intermediate-redshift 
LIRGs (see e.g., Hammer et al. 2005; Marcillac et al. 2006; 
Caputi et al. 2008). Moreover, it has been suggested that LIRGs 
represent phases in the life of galaxies with episodic and ex- 
tremely efficient star formation (Hammer et al. 2005). This sug- 
gests that using one single stellar population (SSP) may not be 
appropriate for modelling individual regions and the integrated 
emission of LIRGs. 



3.1. Modelling of the stellar continuum 

For the modelling of the stellar continuum we used the Bruzual 
& Chariot (2003, BC03 hereafter) models with solar metallicity, 
instantaneous star formation, and a Salpeter (Salpeter 1955) ini- 
tial mass function (IMF) with lower and upper mass cutoffs of 
m\ =0.1 M Q and m u = 100 M Q , respectively. These models have 
a spectral resolution of 3 A across the whole wavelength range 
from 3200 to 9500 A. Using these assumptions we generated 
template spectra covering ages of between 1 Myr and lOGyr. 
The outputs of these models are normalized to a total mass of 
1 M formed in the burst of star formation. The dust attenua- 
tion was modelled using the Calzetti et al. (2000) extinction law, 
which is appropriate for starburst galaxies. 

To fit the optical spectra of the selected regions (see 
Section 2.3) we used the CONFIT code (Robinson et al. 2000), 
which assumes two stellar populations plus a power law in some 
cases. Briefly, CONFIT fits the continuum shape of the ex- 
tracted spectra using a minimum^ 2 technique. For each spec- 
trum, CONFIT measures the flux in ~ 50 - 70 wavelength bins, 
chosen to be as evenly distributed in wavelength as possible, and 
to avoid strong emission lines and atmospheric absorption fea- 
tures (see Rodriguez Zaurin et al. 2009 for details). The relative 
flux calibration error of 10% measured was assumed during the 
modelling. 

For this work we used a large number of combinations of two 
stellar populations to determine which stellar populations dom- 



inate the optical emission of our galaxies. We divided the two 
stellar populations into an ionizing stellar population (< 20 Myr, 
in intervals of 1, 2, ... 10, 20 Myr) with a varying reddening 
(E(B - V) young < 2.0 in increasing steps of 0.1) and an evolved 
stellar population with ages between 100 Myr and lOGyr (100, 
300, 500, 700 Myr; 1, 2, ... 5Gyr and 10 Gyr) with moder- 
ate reddening (E(B - V) evolved = 0, 0.2, 0.4). This preferential 
dust extinction is based on the scenario where the youngest stel- 
lar populations are still partially embedded in their dusty birth 
places, whereas the more evolved stellar populations have al- 
ready moved away from their natal clouds (Calzetti et al. 1994). 
Poggianti & Wu (2000) demonstrated that this scenario was 
compatible with the observed optical spectra of infrared bright 
galaxies. 

The choice of an ionizing population is driven by the fact that 
all the spectra modelled in this paper show Ha in emission with 
EW(Ho') em > 7 A (except for the nuclear region of IC 860, see 
Tables 3 and 4), which sets an upper limit to the age of approxi- 
mately 10-20 Myr for an instantaneous burst (see e.g., Leitherer 
et al. 1999 and below). We did not include very old stellar pop- 
ulations (> 10 Gyr) because these do not appear to have a strong 
contribution to the optical light of local LIRGs and ULIRGs 
(see e.g., Chen et al. 2009; Rodriguez Zaurin et al. 2009). We 
assumed that both stellar populations were formed in instanta- 
neous bursts, which seems an appropriate assumption since the 
physical scales of the selected regions are a few hundred parsecs 
(see Table 3) and contain, at most, a few H n regions, as shown 
in Paper I. 

A priori, fits with x 2 ed < 1 should be considered accept- 
able fits to the overall shape of the continuum (see discussion 
in Tadhunter et al. 2005). However, the absolute value of x 2 ei is 
strongly dependent on the estimated errors. We found that com- 
binations with^- 2 ed > 0.2 produced poor fits to the overall shape 
of the continuum. Moreover, for most regions we found accept- 
able solutions for^ 2 ed < 2 x^^, where x 2 min was me minimum 
value of x 2 ed for a given region. Out of these solutions, we se- 
lected the best fitting models based on a visual inspection of 
the fits to those absorption features with relatively little emis- 
sion line contamination. These include high order Balmer lines, 
the Call K ,13934 line, the G-band /14305, and the Mglb /15173 
band. Finally, we rejected solutions with ages of the young ioniz- 
ing stellar populations older than the upper limits (that is, before 
subtracting the stellar continuum produced by the evolved stars) 
set by the EW(Ha) em values (see Section 3.2). 
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Fig. 5. PMAS spectra of selected regions (see Figure 3 and 
Table 3), arbitrarily scaled, of the three galaxies observed under 
photometric conditions: NGC 23, NGC 2388, and NGC 7771. 
The color used for plotting the spectra corresponds to those of 
Figure 4 (lower panel). 



3.2. Modelling of the hydrogen recombination emission lines 

To model the properties of the hydrogen recombination emis- 
sion lines, we used the Starburst99 model (Leitherer et al. 1999) 
with the same IMF and metallicity assumptions as above to gen- 
erate the time evolution of EW(Ha) em (and also for H/3) for an 
instantaneous burst of star formation. This model is better qual- 
ified for the modelling of populations containing hot massive 
stars (see Vazquez & Leitherer 2005 for details). The EW of the 
Ha emission line resulting from subtracting the modelled con- 
tinuum arising from the non-ionizing stellar population from the 
observed spectra can be used to put further constraints on the 
age (f ne b) of the ionizing stellar populations. Additionally, we 
used the Ha/H/3 emission line ratio measured after subtracting 
the stellar continuum (from both the ionizing and non-ionizing 
stellar populations) to provide an independent estimate of the ex- 
tinction to the gas (E(B-V) ne b). This gas is ionized by the young 
stellar population assumed in the previous section. 



4. The maps of the EW(Ha) em and optical line ratios 

The maps of EW(Ha) em for our sample of LIRGs, covering on 
average the central ~ 4.7 kpc are shown in Figure 1 . The map 
of NGC 7771 covers the central -7.7 kpc X 4.4 kpc. To first or- 
der the EW of the nebular Balmer emission lines can be used as 
indicators of the age of the ionizing stellar populations. The val- 
ues of EW(Ha) em in our sample of LIRGs (see Figure 1), except 
for the nuclear region of IC 860 (Table 4), indicate ages of the 
young stellar populations of between 5 and ~ 10 - 20Myr (see 
e.g., Leitherer et al. 1999, and Section 3.1). For the majority of 
the LIRGs in this sample, as well as for our VLT/VIMOS LIRGs 
of Rodriguez Zaurfn et al. (2010b), the largest values of the EW 
of Ha are not coincident with the peak of the optical continuum 
(the nucleus). However, it is important to note that the EW of the 
Balmer nebular emission lines are also sensitive to the mass of 
the underlying non-ionizing population. In this sense, the EW of 
the Balmer emission lines also provide an estimate of the ratio of 
the current star formation rate compared with the averaged past 
star formation, that is, the burst strength (see e.g., Kennicutt et 
al. 1987; Alonso-Herrero et al. 1996), whether this refers to the 
integrated emission of a galaxy or to individual regions within 
galaxies. In cases of small burst strengths, the observed values 
of the EW only provide upper limits to the age of the current 
star formation burst. Thus, the most likely explanation for the 
smaller nuclear EW(Ha) em , when compared to those of circum- 
nuclear H n regions observed in some galaxies, is a larger con- 
tribution from the underlying (more evolved) stellar population 
(see Kennicutt et al. 1989) and/or a slightly more evolved stellar 
population. 

The regions with the largest EW(Ha) em show values of 
the [Nii]/16584/Ha and [S n]/U67 17,673 1/Ha line ratios (see 
Figure 1) typical of Hn regions in normal star-forming galax- 
ies. In contrast, most nuclear regions in our LIRGs tend to 
show slightly larger lines ratios than the Hn regions of the 
same galaxy. This is clearly seen in the diagnostic diagrams of 
Figure 2 where for each galaxy we plot the spatially resolved (on 
a spaxel by spaxel basis) line ratios as a function of the observed 
value of EW(Ha) em in emission. This is in line with findings for 
the nuclei and H n regions in normal star forming galaxies (see 
e.g., Kennicutt et al. 1989). We note, however, that the line ratios 
of regions with low EW(Ha) em will have the largest corrections 
for the presence of stellar Balmer absorption lines as we shall 
see in Section 6.2. 
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D n (4000) 

Fig. 6. D„(4000) vs. H6a diagram. The small symbols are the 
measurements from the average spectrum of the regions with 
low EW(Ho')e m of each LIRG, whereas the large symbols are 
the measurements corresponding to the nuclear regions (see 
Table 4). For NGC 7591 we plot the values corresponding to 
the integrated spectrum. Note that in the majority of nuclei the 
H6a indices are lower limits due to possible contamination from 
the nebular HS emission line. The black solid line is the time 
evolution (from 1 Myr to 5 Gyr) as predicted by the BC03 mod- 
els using solar metallicity, a Salpeter IMF, and an instantaneous 
burst of star formation. The dashed lines are combinations of dif- 
ferent evolved populations (ages 700 Myr, 1 Gyr, 2 Gyr, 3 Gyr, 
and 5 Gyr) and a young stellar population of 10 Myr. The dotted 
lines represent the fraction in mass of young stars with values of 
0.001, 0.002, 0.005, 0.008, 0.01, 0.02, 0.05, 0.08, 0.1, 0.2, 0.5, 
and 0.8, from right to left. 



The presence of extra-nuclear regions with enhanced val- 
ues of the [Nn]^6584/Hor and [S n]/U6717,6731/Ha ratios rel- 
ative to those of H n regions is again a common property not 
only of the LIRGs studied here, but also of our VLT/VIMOS 
sample of LIRGs (see Monreal-Ibero et al. 2010a) and our 
WHT/1NTEGRAL sample of ULIRGs (see Garcfa-Marfn et al. 
2010). In our sample of LIRGs these regions tend to be as- 
sociated with diffuse emission rather than with high surface 
brightness H n regions, as is also the case for normal and star- 
burst galaxies (Wang, Heckman, & Lehnert 1998). Moreover, 
for our VLT/VIMOS sample of LIRGs Monreal-Ibero et al. 
(2010a) found a correlation between the enhanced optical ra- 
tios and increasing gas velocity dispersion in interacting and 
merger LIRGs, and this correlation is attributed to the presence 
of shocks associated with the interaction processes. Our PMAS 
sample is mostly composed of isolated galaxies and weakly in- 
teracting galaxies (see Paper I), and thus it is unlikely these pro- 
cesses are responsible for the enhanced line ratios. However, 
since the majority of the regions with enhanced [Nn]/16584/Ha 
and [S ii]/t/l6717,6731/Ha' ratios in our sample are observed in 
regions of relatively low EW(Hff) em , we will postpone the dis- 
cussion of this issue after the line ratios are corrected for the 
presence of underlying Balmer stellar absorption features (see 
Section 6.2). 



D n (4000) 



A NGC33 




NGC33B8 
q MCG+03-20-003 
300Myr X IC860 




D n (4000) 



Fig. 7. Same as Figure 6, but showing the effects of reddening 
the young stellar population with E(B - V) young = 0.2 (upper 
panel), E(B - V) young = 0.5 (middle panel), and E(B - V) young 
I (bottom panel), and using the Calzetti et al. (2000) extinction 
law. 

5. The Stellar Populations of local LIRGs 

5.1. Results using the 4000 A break and Balmer line 
absorption features 

Kauffmann et al. (2003a) used a method based on the 4000A- 
break and the Balmer H6 feature in absorption to constrain the 
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mean (light-weighted) age of the stellar population of a galaxy 
and the mass fraction formed in recent bursts of star formation. 
Therefore, this method provides information about the light- 
averaged properties of the stellar populations. From the stel- 
lar continuum spectra generated with BC03, we measured the 
D„(4000) and the U6\ indices (see Section 2.3, and Table 2 for 
definitions of the indices). For an instantaneous burst of star for- 
mation the D„(4000) index increases monotonically with the age 
of the stellar population, while the depth of the H6 line absorp- 
tion feature increases until about 300 - 400 Myr after the burst, 
and then decreases again (black line in Figure 6). Thus the pres- 
ence of high order Balmer lines in absorption is usually inter- 
preted as the signature of an intermediate age (100 Myr - 1 Gyr) 
stellar population. The evolution of these two indices is shown 
in Figure 6 using the BC03 models and is similar to the results 
of Kauffmann et al. (2003a) and Gonzalez-Delgado et al. (2005). 

As discussed in Section 3, we need a combination of (at 
least) two stellar populations to reproduce the observed prop- 
erties of LIRGs, including the D„(4000) and H<5 A indices. This 
is clear from Figure 6, where a single stellar population formed 
in an instantaneous burst (black line) does not reproduce the ob- 
served indices of the nuclei and regions of low-EW(Ha) em in 
our sample of LIRGs. A model with a constant star formation 
rate predicts D„(4000) < 1 .2 for all ages (see figure 21 of Caputi 
et al. 2008), and thus it is not appropriate for our galaxies either. 

Figure 6 shows the result of combining two stellar popula- 
tions. In this diagram, the choice of the age of the ionizing stellar 
population is not critical because D„(4000) and HSa do not vary 
much during approximately the first ~ 10 - 20 Myr of the evo- 
lution of a single stellar population formed in an instantaneous 
burst (see e.g., Gonzalez-Delgado et al. 2005). Thus with this 
kind of diagrams we cannot constrain the age of the youngest 
stellar populations in LIRGs. As we shall see in Section 5.2, 
the combined modelling of the stellar continuum and the neb- 
ular emission lines puts strong constraints on the properties of 
the ionizing stellar population. In Figure 6 we then combined a 
10 Myr population with evolved stellar populations with ages of 
700 Myr, 1 Gyr, 2 Gyr, 3 Gyr, and 5 Gyr. The ages of the evolved 
stellar population are based on the location in this diagram of the 
observed values for the nuclei and regions of low-EW(Ha) em in 
our sample of LIRGs. A scenario where the evolved stellar pop- 
ulation was formed in an instantaneous burst and the current star 
formation is taking place at a constant rate (see Sarzi et al. 2007 
for D„(4000) vs. H£ A diagrams generated under this assumption) 
would underpredict the strength of the H6 absorption feature for 
most of the selected regions in our LIRGs, as was the case for 
the star-forming regions in nuclear rings in the sample of Sarzi 
etal. (2007). 

From Figure 6 it is clear that the main effect of combining a 
young population and an evolved 1-2 Gyr stellar population is 
the change in the observed value of the D„(4000) index, which 
becomes smaller as the fraction in mass of the young stellar pop- 
ulation increases. For a combination with a younger evolved stel- 
lar population (~ 700 Myr) the effect is observed in both indices, 
as is the case for older evolved populations (> 2 Gyr). Based on 
this figure and taking into account that the measured HSa are 
lower limits, the age of the evolved stellar population is between 
1 and 3 Gyr for most of the galaxies in our sample. The only 
exception is NGC 7771 that appears to show the presence of 
an intermediate age (700 Myr- 1 Gyr) stellar population. We also 
show in Figure 6 a range of mass fractions for the young stellar 
population, which in general are relatively small for the nuclei 
and regions of low-EW(Ha) em . 



In Figure 7 we show the effects on the D„(4000) vs. H6a 
diagram from the combination of a reddened young stellar pop- 
ulation and an unreddened evolved population. Figure 7 clearly 
demonstrates that not accounting for the extinction to the young 
stellar population would make us underestimate its mass fraction 
as well as underestimate the age of the evolved stellar popula- 
tion. However, it is also apparent from Figures 6 and 7, that it 
is not possible to disentangle the effects of extinction, ages of 
the stellar populations, and mass contributions from this kind of 
diagrams alone. 

5.2. Results of the modelling of the stellar continuum and 
nebular emission 

As explained in Section 2.3, we selected a number of regions 
in our sample of LIRGs for the study of their stellar popula- 
tions. These include the nuclear regions, the average spectra of 
regions of low-EW(Ho') em , as well as a number of bright H n 
regions in NGC 23, NGC 2388, and NGC 7771. We excluded 
from this analysis three galaxies for various reasons. The nu- 
clear optical spectra of MCG +12-02-001 appears to be com- 
pletely dominated by a young ionizing stellar population (that 
is, we do not see any evidence for the presence of absorption 
features), and thus it was not possible to constrain the evolved 
stellar population. The signal-to-noise ratio of the nuclear spec- 
trum of UGC 1845 did not allow us to constrain the properties 
of the stellar populations. Finally we excluded NGC 7469 be- 
cause of the possible contamination of the optical spectra from 
the AGN non-stellar continuum. We refer the reader to Diaz- 
Santos et al. (2007) for a detailed study of the stellar populations 
of the ring of star formation of NGC 7469 using high angular 
resolution HST photometric data. 

Table 5 summarizes for each galaxy and region the accept- 
able ranges of ages and extinctions for the evolved (f evo ived and 
E(B-V) evo i ved ) and ionizing (f young and E(B-V) young ) stellar pop- 
ulations derived from the stellar continuum fit as explained in 
Section 3.1. We also list in this table the fraction of the light in 
the normalising bin emitted by the ionizing (young) stellar com- 
ponent (/nbX for the models that produce adequate fits. Because 
of the different shapes of the extracted spectra it was not always 
possible to select the same normalizing bin. However, the nor- 
malizing bin was always selected to be located within the wave- 
length range of 4400 - 4800 A and usually spanning ~ 100A. 
Figure 8 presents a few examples of fits to the stellar continuum 
for a number of regions in our sample. 

For the majority of the regions studied here except for the 
optical nucleus of NGC 7771 (see Section 5.2.1), we find that 
stellar populations with ages between 100 Myr and 500 Myr do 
not make a strong contribution to the optical light. Ages for 
the evolved stellar populations of between 0.7 and 5 Gyr (or 
even 10 Gyr in some cases) provided reasonable fits to the op- 
tical continua. These ages are within the range of ages derived 
for massive spiral galaxies (Gallazzi et al. 2005 and references 
therein). Another noteworthy result is that the regions with low- 
EW(Ha) em in the galaxies tend to have slightly older ionizing 
stellar populations and less extinction than other regions for the 
same galaxy. The ages of the evolved stellar populations agree 
with our findings in the previous section using the D„(4000) and 
H<5a diagram, keeping in mind that H£ A may only provide a 
lower limit to the age of the evolved stars for some of the re- 
gions. We finally note that for models including evolved plus 
young stellar components it is particularly hard to distinguish 
between ages in the range 2 - 10 Gyr. Therefore, in some cases 
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Fig. 8. (a) Examples of fits to the stellar continuum of regions in our sample of galaxies. The bluest rest-frame wavelength used 
for the modelling is 3800 A. For each region, the left panel shows the observed spectrum (thick color line) and the model spectrum 
(thin black line) in arbitrary units. The spectra are normalized to unity at a wavelength within the normalizing bin (4400 - 4800 A). 
The model parameters for each region are given at the top right of each plot. The right panels are blow-ups of some spectral regions 
of interest. We mark the high order Balmer lines as well as HS, Hy, and, Hfi as dotted lines. The Can H and Can K lines (lower 
right panel), the G-band (middle right panel), and the Mgib band (upper right panel) are marked as thick solid lines. The spectral 
resolution of the BC03 models has been slightly degraded to match approximately that of our PMAS spectra. 
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Table 5. Ranges of ages, extinctions and light fractions from young stars as derived from the stellar continuum modelling. 



Galaxy 


Region 


^evolved 

Gyr 


E (B - V )evol ved 


Myr 


E(B- V) vn i ln ,r 


,/nb 

% 


NGC 23 


Nucleus 


0.7-5 


0-0.4 


< 10 


0.4-1.2 


15-56 




Low-EW(Ha) em 


0.7- 1 


0-0.4 


< 20 


0.4-1.0 


20-55 




HII-1 


2-10 


0-0.4 


< 6 


0.4-0.9 


30-60 




T4TT 9 


fl 7 - 1 
u. / 1 


o ? - n 4 


A — 7 




90 — 49 




HII-3 


0.7-1 


0.4 


<7 


0.4 


44-57 




HII-4 


\J. J 1 




v o 


9 - f> 








9 - s 


4 


(S - 7 


9—10 






Low-EW(Ha) em 


1 


0.4 


7-20 


0.7 - 0.9 


26-40 




HII-1 


0.7- 1 


0-0.4 


6-7 


0.8-1.0 


63-74 




HII-2 




4 


u 


1 3 


40 


MCG +02-20-003 


Nucleus 


0.7-2 


0.2-0.4 


< 6 


1.5 - 1.8 


43-73 


IC 860 


Nucleus 


0.5-1 


0.2-0.4 


< 20 


0.3-1.4 


18-53 


NGC 5936 


Nucleus 


1-5 


0-0.2 


5-7 


1.0-1.1 


76-84 




Low-EW(Ha) em 


0.5 - 10 


0-0.4 


7-10 


0.3-0.8 


52-71 


NGC 6701 


Nucleus 


0.7-5 


0-0.4 


< 7 


0.4-0.5 


43-68 




Low-EW(Ha) cm 


0.7-2 


0-0.2 


7-10 


0.2-0.4 


52-71 


NGC 7591 


Nucleus 


2-10 


0-0.4 


< 7 


0.4-1.5 


20-63 


NGC 7771 


Nucleus 


0.3 - 0.5 


0.2 


4-7 


0.3-0.6 


37-53 




Nucleus2 


2-5 


0.4 


8-20 


0.4-0.6 


48-63 




Low-EW(Ha) em 


0.5 - 0.7 


0.4 


6-10 


0.4-0.6 


45-57 




HII-1 


0.7 - 10 


0-0.4 


< 6 


0.4-0.9 


50-75 




HII-2 


2-10 


0-0.4 


< 7 


0.4-1.1 


30-70 




HII-3 


2-10 


0-0.4 


< 8 


0.3-1.1 


25-58 



Table 6. Examples of stellar continuum and nebular modelling. 



Galaxy 


Region 




Stellar continuum parameters 




Nebular fit 






^evolved 


E(B - V) evo i ve( j 


^young 


E(B-V) young 


/nb 


^neb 


E(B-V) neb 






Gyr 




Myr 




% 


Myr 




NGC 23 


Nucleus 


2 





7 


0.7 


31 


6.9 


0.6 


NGC 23 


HII-4 


1 


0.2 


5 


0.2 


44 


5.7 


0.5 


NGC 2388 


HII-1 


1 





6 


1.0 


80 


6.0 


1.0 


IC 860 


Nucleus 


1 


0.4 


9 


0.3 


27 


8.3 




NGC 6701 


Low-EW(Ha) em 


2 





8 


0.2 


58 


7.6 


0.3 


NGC 7771 


Nucleus 


0.5 


0.2 


7 


0.3 


40 


6.1 


0.5 



we are not able to put strong constraints on the age of the stellar 
populations for ages older than ~ 1 - 2 Gyr. 

As can be seen from Table 5 the ionizing stellar populations 
always contribute a minimum of / NB > 15% to the optical emis- 
sion in all the modelled spectra, and they dominate the emission 
(/nb ^ 50%) in, at least, 7 of the 23 spectra modeled. On the 
other hand, in the nuclear regions of NGC 23 and IC 860, and 
NGC 777 1 , the evolved stellar populations may be the main con- 
tributors to the optical light (see also Table 6). Our spectra do 
not sample the near-UV spectral region, which is important for 
constraining the properties of the ionizing stellar populations. 
Therefore, the ages of such stellar populations are not well con- 
strained when fitting the stellar continuum alone. 

We were able to put tighter constraints on the young stel- 
lar populations when we included the nebular fitting. Table 6 
gives the parameters of the stellar continuum and nebular fits 
for those regions shown in Figure 8. The nebular ages and 



extinctions (f ne b and E(B-V) ne b) were derived as explained in 
Section 3.2, for an acceptable combination of stellar populations 
close or at the minimum value of ^ .. We find ages of the ion- 
izing stellar populations from the EW of Ha of between 5.6 and 
8.8 Myr. The extinctions to the ionizing stellar stars range be- 
tween E(B - V) = 0.2 and E(B - V) = 1.8 (Tables 5 and 6). 
The extinctions to the ionizing stellar populations are always sig- 
nificantly greater than those to the evolved stars, and generally 
consistent with those derived to the young stars from the stellar 
continuum modelling. 

Our results are, in a broad brush sense, consistent with those 
of Rodriguez Zaurin et al. (2009) for their sample of ULIRGs. 
That is, the optical spectra can be modeled using a combina- 
tion of an evolved plus a young stellar population. However, the 
extracted spectra of our LIRGs show, in most cases, deeper G 
and Mgib bands than those of the ULIRGs of the Rodriguez 
Zaurin et al. (2009) sample. This suggests that the evolved stel- 
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Fig. 8. (b). Same as Figure 8a. For the nuclear region of NGC 7771 we also show the positions of the Hei features at 3820, 4026, 
4388, and 4471 A. 



lar populations are somewhat older for the LIRGs in our sam- light of ULIRGs could be the result of enhanced star forma- 

ple. In fact, Rodriguez Zaurin et al. (2009) found adequate fits tion coinciding with the first pass of the merging nuclei, along 

for their ULIRGs which included evolved stellar populations with a further, more intense, episode of star formation occur- 

of 0.3 - 0.5 Gyr, while this was rarely the case for our sam- ring as the nuclei finally merge together (Rodriguez Zaurin et 

pie of LIRGs. The stellar populations dominating the optical al. 2010a). Our sample of LIRGs on the other hand, is mostly 
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composed of relatively isolated spiral galaxies and weakly inter- 
acting galaxies (see Paper I), with moderate IR luminosities (see 
Section 2.1). Given the relatively small fraction of strongly inter- 
acting/merger systems in our sample compared to the Rodriguez 
Zaurfn et al. (2009) ULIRG sample, it may not be unexpected 
that intermediate-age (~ 100 - 500 Myr) stellar populations do 
not dominate the optical light of our LIRGs. The only exception 
in our sample are the central regions of NGC 7771 (see discus- 
sion in Section 5.2.2). Moderate luminosity spiral-like LIRGs 
may be constantly forming stars and may have not undergone a 
major burst of star formation in the last 1-2 Gyr, as is the case of 
normal spiral galaxies (Kauffmann et al. 2003b). Our results are 
also in accord with the observational findings of Poggianti & Wu 
(2000) that most isolated systems in their sample of IR-bright 
galaxies showed on average more moderate Balmer absorption 
features than the interacting systems. The effects of possible mi- 
nor mergers are likely to be difficult to evaluate using the stel- 
lar populations as models and observations suggest that it is the 
satellite galaxy rather than the primary galaxy that is more sus- 
ceptible to enhanced star formation (Woods & Geller 2007; Cox 
et al. 2008). 

5.2.1 . The central regions of NGC 7771 

There is dynamical evidence that NGC 7771 is weakly inter- 
acting with NGC 7770 (Keel 1993) and is located in a group 
of galaxies. It is then possible that the interaction process re- 
sulted in a strong burst of star formation in the past. For in- 
stance, the ring in this galaxy, clearly detected in our Paa images 
(see Paper I), appears to have a complex star formation history 
with evidence for multiple generations of stars (Davies, Alonso- 
Herrero, & Ward 1997; Smith et al. 1999; Reunanen et al. 2000). 

The optical nucleus of NGC 777 1 is the only region in our 
sample of LIRGs for which adequate fits were obtained with a 
large contribution to the optical light from an intermediate-age 
(300 - 500 Myr) stellar population (see Table 5 and Figure 8). 
The fit to the average spectrum of regions with low EW(Ha) em 
also required an evolved stellar population of 500 - 700 Myr. 
In the case of the nuclear spectrum there is evidence for an 
additional stellar population based on the presence of He i ab- 
sorption features at various wavelengths, which are indicated in 
Figure 8.b. These absorption lines are strongest for stellar pop- 
ulations with ages in the range 20 - 50 Myr (Gonzalez-Delgado 
et al. 1999; 2005), and are not observed in stellar populations 
older than 100 Myr (the lifetime of B stars). The presence of this 
important population of non-ionizing stars dominated by B stars 
was already infered by Davies et al. (1997). Given the clear ev- 
idence for the presence of B stars, for this region we also tried 
combinations including stellar populations of 30 - 80 Myr for the 
young component. We found that a combination with a dominant 
(/nb = 70 - 80%), low reddening (E(B - V) young ~ 0.2) stellar 
population of 40 - 50 Myr plus an unreddened evolved stellar 
population of few Gyr provided an acceptable fit. We note that 
with this combination the stellar population responsible for ion- 
izing the gas is not accounted for. Therefore, in this particular 
case, it is likely that models including a larger number of stellar 
components (at least three) would be more adequate. We again 
emphasize that the nuclear region of NGC 777 1 is the only re- 
gion where we find a clear evidence for these He i features. In 
this respect, it is not clear if the optical nucleus of this galaxy 
may be a special region in this galaxy (see Davies et al. 1997 for 
a discussion on this issue). It is important to recall at this point 
that the optical nucleus of NGC 777 1 is probably not the true 
nucleus of the galaxy, as it does not coincide with the center of 
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D„(4000) vs. Ha a (lower panel) diagrams. The black solid line 
is the time evolution (from 1 Myr to 5 Gyr) as predicted by the 
BC03 models using solar metallicity, a Salpeter IMF, and an in- 
stantaneous burst of star formation. The dashed lines are com- 
binations of different evolved populations (ages 700 Myr, 1 Gyr, 
2 Gyr, 3 Gyr, and 5 Gyr) and a young ionizing stellar population 
of 6 Myr. The dotted lines represent the contribution in mass of 
the two stellar populations, as in Figure 6. 

the bright ring of star formation, the peak of near-IR emission 
(see Paper I), or even with the region with the largest value of 
EW(Ha) em (see Figure l.j). 

5.3. Predictions for the H/3 and Ha stellar absorption features 

To study in detail the optical line ratios and in particular those 
of the diffuse regions, we need to correct for the presence of H/3 
and Ha stellar absorption features. However, given the strong 
star formation activity in our sample, it is difficult to measure 
reliably these absorption features because of the strong contami- 
nation produced by the Balmer recombination lines in emission. 
An alternative approach is to generate theoretical diagrams of 
D„(4000) vs. Ha A and D„(4000) vs. H/3 A measured from the 
spectra generated from the combination of the BC03 young and 
evolved stellar populations as explained in Section 3.1. The H/3 a 
and HffA indices for the stellar absorption features are defined 
in a similar way to that used for the H6 absorption feature (see 
Section 2.3), and thus have a positive value. The line windows, 
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and the blue and red pseudo-continuum windows defined by 
Gonzalez Delgado et al. (2005) are given in Table 1. For the 
ages of the evolved stellar population we used those derived 
in Section 5.2. For the ionizing stellar population we chose an 
age of 6 Myr, which is representative of our sample of LIRGs 
(Tables 5 and 6), although the results are not strongly dependent 
of chosen age of this population. For the mass fractions we used 
the same values as those shown in Figure 6. 

From the D„(4000) vs. H/3a diagram (Figure 9, upper panel), 
we can see that for ages of the evolved stellar population of be- 
tween 1 and 5 Gyr, the predicted average value of the H/3a index 
is ~ 5±1 A. This value is almost independent of the mass fraction 
in young stars, except for cases where the stellar mass is dom- 
inated by the contribution from ionizing stars. For NGC 7771, 
which is the clearest case in our sample for the presence of an 
intermediate age stellar population, the predicted value would be 
HjSa ~ 7 + 1 A. In this case, the predicted value is more sensi- 
tive to the mass in young stars and the age of the evolved stellar 
population. For comparison, the measurements of H/3a reported 
by Kim et al. (1995) for the central 2 kpc of the galaxies in com- 
mon with our sample are 3 and 6 A, although Kim et al. (1995) 
pointed out these values were admittedly subjective because of 
the method they used for fitting their data. 

Figure 9 (lower panel) shows the predictions for the 
D„(4000) vs. Hq-a diagram. In the case of Hcka there is a very 
small dependence of the predicted value with the mass fraction 
in young stars. For the majority of our LIRGs with ages of the 
stellar population 1-5 Gyr the average predicted value for the 
index would be HffA ~ 3 + 0.5 A, whereas for NGC 7771 we 
would predict Ha a ~ 4 + 0.5 A. For comparison, Moustakas & 
Kennicutt (2006b) found an average Ha stellar absorption cor- 
rection of 2.8 + 0.4 A for the integrated spectra of a sample of 
nearby star forming galaxies. 



6. Excitation Conditions in local LIRGs 

6.1. Spatially resolved diagnostic diagrams 

Diagnostic diagrams using bright optical emission lines (Balwin 
et al. 198 1 and Veilleux & Osterbrock 1987) are useful for differ- 
entiating between the various sources of excitation of the gas in 
the nuclei of galaxies and their integrated emission. With optical 
IFS data we can additionally study the distribution of the ioniza- 
tion structure of spatially resolved regions (that is, on a spaxel- 
by-spaxel basis) of nearby galaxies (e.g., Garcia-Mann et al. 
2006; Garcfa-Lorenzo et al. 2008; Blanc et al. 2009; Stoklasova 
et al. 2009; Monreal-Ibero et al. 2010b; Garcia-Mann et al. 
2010), as well as of Hn regions in nearby galaxies (e.g., Relano 
et al 2010). The diagnostic diagrams for the spatially resolved 
measurements for each of the LIRGs in our sample are shown in 
Figure 2. The measurements for the individual spaxels are color 
coded according to the observed value of the EW(Ha) em in emis- 
sion. In these diagrams we plotted the empirical and theoretical 
boundaries derived by Kauffmann et al. (2003c) and Kewley et 
al. (2001). These boundaries are shown for reference as they may 
provide clues about the dominant excitation mechanism: ioniza- 
tion by young stars, shocks, or AGN ionization. 

The first result worth noticing is that a large number of 
the spatially resolved measurements in the [Om]/15007/HjS vs. 
[Nn]/165 84/Ha diagram (Fig. 2, left panels) fall in the com- 
posite region, and in particular a large fraction of those spaxels 
with EW(Ha) em < 20 A. A similar result is found for ULIRGs 
(Garcia-Mann et al. 2010). On the other hand, most of the spax- 



els with EW(Ho-)em > 60 A are located in the H n region of the 
diagram, with the only exception of some regions in NGC 7469. 
There are also differences from galaxy to galaxy. For instance 
the majority of the spaxels of MCG +12-02-001 and NGC 5936 
are located in the Hn region of the diagram, while other cases 
such as UGC 1845 most spaxels are in the composite region. 
The composite region on this diagram lies between the observa- 
tional AGN/H ii boundary and the "maximum starburst line" of 
Kewley et al. (2001), above which the observed line ratios can- 
not be explained by star formation alone. Kewley et al. (2006) 
interpreted the observed ratios in the composite region as pro- 
duced by a metal-rich stellar population and an AGN (but see 
Cid Fernandes et al. 2010 and references therein, for an oppos- 
ing point of view). While the Kewley et al. argument is valid for 
nuclear and even integrated spectra of galaxies, in our sample a 
large fraction of spaxels with enhanced line ratios are detected 
in the extra-nuclear regions of galaxies without an AGN. 

The number of spatially resolved measurements located in 
the LINER region of the diagrams of the other two diagnostic 
diagrams (Fig. 2, middle and right panels), is smaller than in 
the diagram with the [Nn]/16584/Ha ratio (see also Veilleux et 
al. 1995). This is in part a sensitivity issue, especially for the 
relatively faint [Oi]/I6300 line, which is, on the other hand, a 
very good shock tracer. Still, a large fraction of spaxels with 
EW(Ha) em < 20 A, mostly in NGC 23, MCG+02-20-003, and 
NGC 7771, are in the LINER region of these two diagnostic 
diagrams. The presence of extra-nuclear regions in LIRGs and 
ULIRGs with LINER-like excitation has often been interpreted 
as an indication for the presence of large scale shocks where 
interactions are playing a major role (Monreal-Ibero, Arribas, 
& Colina 2006; Monreal-Ibero et al. 2010a) or to shocks due 
to the presence of outflowing nuclear gas (see e.g., Armus et 
al. 1989). Calzetti et al. (2004) using a [Om]A50Q7/H/3 vs. 
[S n]/L167 17,673 l/Ha diagram with spatially resolved measure- 
ments of four nearby starburst galaxies, found evidence for the 
presence of non-photoionized gas. These authors also demon- 
strated that shocks from supernovae and stellar winds are able 
to provide sufficient mechanical energy to account for the non- 
photoionized gas in their galaxies. In the case of this sample of 
LIRGs, most galaxies do not show evidence for strong interac- 
tions either from their morphologies or from the observed Ha ve- 
locity fields (see Paper I). Our LIRGs, however present very high 
central Ha surface brightnesses and thus star formation rates 
per surface area (see Alonso-Herrero et al. 2006), which would 
imply a contribution from supernovae. One important caveat to 
keep in mind is that the observed line ratios in the diagnostic 
diagrams of Figure 2 have not been corrected for the presence 
of stellar absorption. These corrections will be most relevant for 
those spaxels with the lowest equivalent widths of the hydrogen 
recombination emission lines, as we shall see in the following 
section. 

6.2. Spatially resolved [Nn]A6584 /Ha vs. 
[Sn]AA671 7,6731/Ha diagrams 

In our own Galaxy and other galaxies enhanced [Nii]/16584/Hq' 
and [S ii]/U6717,6731/Hq' line ratios appear to be associated 
with the presence of diffuse ionized gas (DIGQ, see the re- 
cent review by Haffner et al. 2009 and references therein). The 
presence of DIG emission has been detected spectroscopically 
in the extra-planar emission of edge-on galaxies (Rand 1996; 



3 This gas is also referred to as diffuse ionized medium or DIM, and 
warm ionized medium or WIM. 
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Fig. 10. Spatially resolved [Nn]/16584/Ha vs. [Sn],Ll6717,6731/Ha diagrams for the full PMAS sample. The individual measure- 
ments (a total of 1852 spaxels) in both panels are color coded according to the observed EW(Ha) em in emission as in Figure 2. The 
left panel shows the ratios not corrected for underlying Ha stellar absorption, whereas the right panel are the line ratios corrected 
for underlying Ha stellar absorption as described in Section 5.3. The contours are plotted to help assess the data point density in 
bins of 0. 1 dex in both emission line ratios. 



Miller & Veilleux 2003). A similar result was infered based on 
the enhanced low-ionization emission relative to that of H n re- 
gions of the integrated emission of galaxies (see e.g., Lehnert & 
Heckman 1994; Wang, Heckman, & Lehnert 1997; Moustakas 
& Kennicutt 2006a) as well as on spatially resolved measure- 
ments (see e.g., Calzetti et al. 2004). In Paper I we already dis- 
cussed the possibility of the presence of this diffuse emission in 
our sample of LIRGs, as in general the integrated line ratios of 
the galaxies are greater than the typical values observed in disk 
Hn regions (see Kennicutt et al. 1989). Moreover, the regions 
of enhanced line ratios in our LIRGs are generally not associ- 
ated with regions of high Paa surface brightness, that is, bright 
H ii regions. Clear examples of this are NGC 23, NGC 7591 and 
NGC 7771 (see Figure 1 and Paper I). Regions of enhanced low 
ionization emission can also be associated to the presence of an 
AGN in the form of an ionization cone, as is the case of one of 
the nuclei of the interacting LIRG Arp 299 (Garcfa-Marin et al. 
2006). 

To study in more detail the excitation conditions in 
our sample of LIRGs, we produced spatially resolved 
[Sn]^6717,6731/Ha vs. [Nn]A6584/Ha diagrams for each of 
the galaxies (see Figure 3). The advantages of these diagrams 
over the standard diagnostic diagrams are twofold. First they 
contain more data points than the diagnostic diagrams (see the 
statistics in Table 1 for the number of spaxels with measurements 
for each line ratio) because Hy8 in LIRGs is heavily affected 
by extinction and stellar absorption. Second, as we showed in 
Section 5.3, the corrections for the presence of stellar absorption 
in Ha are less dependent on the results of the stellar population 
models than those for H/3. 

A comparison between the observed [S n]AA61 17,673 1/Ha 
vs. [Nn]/16584/Ha ratios and predictions from different models 
(ionization by young stars, shocks, AGN photoionization) can 
shed some light on the dominant excitation conditions in our 
sample of LIRGs. In Figure 3 we show the Dopita et al. (2006) 
models for evolving Hn regions. In these models the ionization 
parameter is replaced by the % parameter, which is defined as 
the ratio of the mass of the ionizing cluster to the pressure of the 



interstellar medium. We chose models with solar and twice so- 
lar metallicity based on the derived abundances of our galaxies 
from the integrated line ratios over the central few kpc (Paper I) 
corrected for stellar absorption (see Section 5.3). To estimate the 
abundances, we used the Pettini & Pagel (2004) empirical cali- 
bration based on the [Om]/i5007/HjS and [Nn]/16584/Ha ratios, 
also known as the 03N2 index (see Alloin et al. 1979). All the 
LIRGs in our sample have near solar or super-solar abundances 
(see Table 7), for a solar abundance of 12 + log(0/H) = 8.66 
(Asplund et al. 2004), and are within the derived abundances of 
the large sample of LIRGs studied by Rupke et al. (2008). 

Although we do not intend to use the Dopita et al. mod- 
els for dating the Hn regions, it is clear that line ratios of 
[S n] /1/167 17,673 1/Ha -0.6-1 could be produced by evolved 
and metal rich Hn regions (see details in Dopita et al. 2006). 
As can be seen from Figure 3, the [S n]AA67 17,673 1/Ha vs. 
[Nn]/16584/Ha diagrams of the central regions of some galax- 
ies (e.g., MCG+12-02-001 and NGC 5936) could be mostly ex- 
plained as emission coming from H n regions, although it has 
to be noted that there are no spaxels with EW(Ha) em < 5 A in 
the central regions of these galaxies. On the other hand, galax- 
ies like NGC 23, NGC 7591 and NGC 7771 show a signifi- 
cant numbers of spaxels with EW(Ha) em < 20 A whose line 
ratios (not corrected for stellar absorption) cannot be explained 
by the Dopita et al. (2006) Hn region models. Shock models 
such as those of Allen et al. (2008) could explain line ratios 
log ([Nn]^6584/Ha) > -0.7 and log ([S n]AA61 17,673 1/Ha) > 
-0.7. These models are not plotted in Figure 3, but see figure 8 
of Monreal-Ibero et al. (2010a) and a discussion by Miller & 
Veilleux (2003). 

To get a more global picture of the excitation conditions in 
LIRGs, in Figure 10 (left panel) we plot a [Nn]/16584/Ha vs. 
[S n] A/167 17,673 1/Ha diagram for the full PMAS sample, show- 
ing a total number of 1852 spaxels again color coded in terms 
of EW(Ha) em in emission. It is clear from this figure that those 
spaxels with the smallest EW(Ha) em tend to show the largest val- 
ues of the line ratios, indicating that in part the enhanced values 
are due to the presence of Ha stellar absorption. The right panel 
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of Figure 10 shows the same diagram after correcting statisti- 
cally each spaxel for the presence of Ha in absorption, with the 
average values of Ha A derived in Section 5.3. The effect of cor- 
recting the emission line ratios for stellar absorption in our sam- 
ple is evident, especially for spaxels with EW(Ha) em < 20 A. 
However, as can be seen from the histograms in Figure 11, even 
after the correction for stellar absorption there is still a signifi- 
cant number of spaxels with enhanced line ratios. For example, 
after correction for stellar absorption approximately 25% of the 
spaxels show log ([S n]AA67 17,673 1/Ha) > -0.2, which is the 
largest line ratio allowed by the 2Z Dopita et al. (2006) models. 
This suggests that there are other mechanisms apart from young 
stars ionizing the gas in our sample of LIRGs. 

We can now compare our observations with those of 
Monreal-Ibero et al. (2010a) for the VLT/VIMOS sample of 
LIRGs. The VIMOS sample was observed with a mean spatial 
sampling of 270 pc per spaxel, which is similar to that of the 
PMAS data. The Monreal-Ibero et al. (2010a) sample contains 
galaxies chosen to probe the different morphologies observed 
in LIRGs from relatively isolated galaxies, to interacting galax- 
ies to mergers. In that study we found that the median value of 
the [Nn]/16584/Ha ratios measured on a spaxel-by-spaxel ba- 
sis has a very weak dependence with the morphology of the 
system. The distributions of the [S n]AA67 17,673 1/Ha and the 
[Oi]/I6300/Ha ratios, on the other hand, tend to be more en- 
hanced for galaxies classified as interacting and mergers. Our 
PMAS sample is smaller than the VIMOS sample, but it is rep- 
resentative of the LIRG class, in the sense that it comprises the 
majority of the IRAS RBGS LIRGs at d < 75 Mpc that can be 
observed from Calar Alto. As such, the PMAS sample is mostly 
composed of LIRGs with relatively low IR luminosities, and thus 
dominated by isolated spiral galaxies or systems in weak inter- 
action (see also Sanders & Ishida 2004). The median values of 
the line ratios (on a spaxel-by-spaxel basis) in the PMAS sam- 
ple corrected for stellar absorption are [Nn]/16584/Ha = 0.45 
and [S n] AA61 17,673 1/Ha = 0.44 (see also Figure 11). For the 
statistics we used only spaxels with measurements of both line 
ratios. These ratios are comparable to, although slightly greater 
than, those measured by Monreal-Ibero et al. (2010a) for the 
VLT/VIMOS sample of LIRGs. The slight difference may be 
due to the fact that the PMAS data probe more regions (spax- 
els) with low EW(Ha) em than in the VLT/VIMOS sample (see 
Rodriguez Zaurin et al. 2010b). This in turn results in the PMAS 
data being more sensitive to diffuse emission. 



6.3. Diffuse emission 

Regardless of the main mechanism responsible for exciting the 
diffuse gas in galaxies (e.g., massive stars, escaping photons 
from H ii regions, shocks), it is clear that quantifying the fraction 
of diffuse emission in galaxies is important for understanding is- 
sues such as the early reionization of the universe. Traditionally 
the fraction of diffuse emission in galaxies was obtained by iden- 
tifying and doing photometry of Hn regions, to separate the 
emission coming from Hn regions from that of the DIG (see 
e.g., Zurita, Rozas, & Beckman 2000; Oey et al. 2007, see re- 
view by Haffner et al. 2009). The main results from this method 
were that the mean fraction of diffuse emission in galaxies was 
50 - 60% with no correlation with the Hubble morphological 
type (Thilker et al. 2002; Oey et al. 2007). However, Oey et al. 
(2007) found that starburst galaxies showed lower fractions of 
diffuse Ha emission when compared to other galaxies. 
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Fig. 11. Histograms showing the distributions of the spatially re- 
solved [Nn]/i6584/Ha (upper panel) and [S n]AA67 17,673 1/Ha 
(lower panel) line ratios for the full PMAS sample, for a total of 
1852 spaxels. The thick line empty histograms are the observed 
line ratios, while the thin line filled histograms are the line ratios 
corrected for stellar absorption. 

Table 7. Abundances and fraction of diffuse ionized emission in 
the central regions of LIRGs. 

DIG fraction 
NICMOS 2xMW - MW 

30 +1 7 3 
2l +9 . 
50+ 8 f3 

S« 

^-10 

76+ 19 
44+1^ 

... -m 

Notes. — The abundances are derived using the 03N2 index measured 

in the integrated PMAS spectra (Section 6.2). The NICMOS DIG frac- 
tions are from H n region photometry on the WST/NICMOS Paa images. 
The 2 x MW and MW DIG fractions are from comparing the spatially 
resolved [S n]/16717/Ha ratios in our LIRGs with those of the Milky 
Way H ii regions and DIG (see Section 6.3 for details). 
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In Alonso-Herrero et al. (2006) we measured the Paa 
fluxes of H ii regions, and estimated the total Paa fluxes using 
HSr/NICMOS narrow band images, for similar FoVs to those 
of the PMAS data. Due to the relatively small FoV of the Paa 
NICMOS images (NIC2 camera 19" x 19") the main uncertainty 
in measuring Paa fluxes, both in H n regions and the total (H n 
+ DIG) emission is the background removal. On the other hand, 
narrow-band Paa imaging does not suffer from problems asso- 
ciated with contamination by the [N n] doublet associated with 
using narrow-band Ha imaging (see Blanc et al. 2009 for a dis- 
cussion on the issue). For each galaxy we measured the total 
Paa luminosity over the NICMOS FoV, and converted it to Ha 
flux assuming case B recombination (Ha/Paa = 8.6, Hummer 
& Storey 1987). Using the total Ha luminosity and the Ha lu- 
minosity in H ii regions given in table 3 of Alonso-Herrero et 
al. (2006), we computed the DIG fraction for each galaxy. The 
DIG fractions are between 8 to 76% (Table 7). The range given 
for each galaxy takes into account the uncertainties in calculat- 
ing the total Paa emission from the HST/NICMOS narrow-band 
imaging (see Alonso-Herrero et al. 2006 for details). The DIG 
fractions given in this table for each LIRG correspond only to the 
central few kiloparsecs, as neither the NICMOS nor the PMAS 
observations cover the full extent of the galaxies. 

Alternatively Blanc et al. (2009) computed the DIG frac- 
tion in M51 from optical IFS by comparing the observed 
[S ii]/l6717/Ha line ratios with the typical values measured for 
H ii regions and the DIG in the Milky Way (MW) by Madsen 
et al. (2006), and taking into account the difference in metal- 
licity between M51 and the MW. The PMAS 1" spaxel covers 
approximately 300 pc for the typical distances of our galaxies 
(Section 2.1). Defining the sizes of Hn regions, even at HST 
resolutions, is not straightforward, but results for LIRGs at dis- 
tances similar to those of the PMAS sample, indicate typical 
sizes of ~ 150 - 200 pc (Alonso-Herrero et al. 2002). We thus 
expect that each measurement of a given line ratio contains both 
H ii region and DIG emission, as assumed by the Blanc et al. 
(2009) method. In the MW the average [S n]i6717/Ha line ra- 
tios of H ii regions and the diffuse emission are 0. 1 1 +0.03 and 
0.34 ±0.13, respectively (Madsen et al. 2006; Blanc et al. 2009). 
In Table 7 we list the DIG fractions we obtained with the line 
[S n]/l6717/Ha ratios corrected for stellar absorption and the 
Blanc et al. (2009) method. For each galaxy we give a range, 
with the smallest value corresponding to scaling the MW Hn 
and DIG ratios to ~ 1.7 times solar, and the largest ratio for so- 
lar metallicity. 

Table 7 shows that the fraction of diffuse emission in the 
central few kpc of our LIRGs varies significantly from galaxy 
to galaxy. Oey et al. (2007) found that starburst galaxies, de- 
fined as galaxies having Ha surface brightnesses above 2.5 x 
10 39 erg s _1 kpc~ 2 within the Ha half radius, have relatively 
small DIG fractions (< 60%), when compared to non-starburst 
galaxies. According to this definition all our galaxies are star- 
bursts (see Alonso-Herrero et al. 2006), and thus in most cases 
our DIG fractions are below 60%. The large DIG fraction in- 
ferred for NGC 7591 may be explained if some of enhanced 
line ratios are produced in a narrow line region associated with 
the AGN in this galaxy. We also find that the fractions of DIG 
emission in the central few kpc obtained with the two methods 
described above mostly agree with each other, within the uncer- 
tainties (Table 7). We note however, that the line ratio method 
always tends to provide a smaller DIG fraction than the H ii re- 
gion photometry method (see also Blanc et al. 2009 for M51). 
The two LIRGs with the most discrepant fractions using the 
two methods are NGC 5936 and NGC 7771. In the case of 



NGC 5936 we attribute the discrepancy to the fact that, except 
for the very nuclear regions, the other H ii regions in the FoV of 
the NICMOS images appear quite diffuse (see Alonso-Herrero 
et al. 2006, and Paper I), and thus the DIG fraction estimated 
from the NICMOS Paa images is clearly just an upper limit. In 
the case of NGC 7771, even if we applied the smaller correc- 
tion for the Ha absorption feature (that is, if the evolved stellar 
population had an age of a few Gyr), the DIG fraction from the 
[S n]/i6717/Ha line ratios would be in the range of 3 to 13%. 
The most likely explanation is that the faintest H ii regions just 
outside the bright ring of star formation were not identified in 
the NICMOS Paa images, and this resulted in an overestimate 
of the diffuse emission. 



7. Summary 

This is the second paper of a series presenting a PMAS IFS 
study of 11 local (average distance of 61Mpc) LIRGs cover- 
ing the central few kiloparsecs of the galaxies. We selected the 
galaxies from the complete volume-limited sample of LIRG of 
Alonso-Herrero et al. (2006), which has an average IR luminos- 
ity of log(Li R /L ) = 11.32. As such, our sample of LIRGs is 
mostly composed of isolated spiral galaxies, galaxies with com- 
panions, and weakly interacting galaxies. The PMAS observa- 
tions covered a spectral range of 3700 - 7200 A, with a spectral 
resolution of 6.8 A (FWHM). We used a l"-size spaxel that pro- 
vides a typical physical sampling of approximately 300 pc. The 
PMAS observations were complemented with our own existing 
ffi>r/NICMOS observations of the near-infrared continuum and 
the Paa emission line. We combined the spectral information 
provided by the PMAS IFS and the spatial information provided 
by the HST/NICMOS observations to study in detail the stel- 
lar populations, the excitation conditions of the gas, and diffuse 
emission in local LIRGs. 

To study the stellar populations we selected the nuclear re- 
gions of each galaxy, as well as a number of bright H n regions 
in NGC 23, NGC 2388, and NGC 7771. We also extracted for 
each galaxy, when possible, the average spectrum of regions 
with low values of EW(Ha) em . These regions were selected such 
that they did not contain high surface brightness H n regions to 
minimize the contamination by hydrogen nebular emission in 
the stellar absorption features. We used BC03 and Starburst99 to 
fit the stellar continuum and the nebular emission, respectively. 
For the modelling we assumed two stellar populations formed 
in instantaneous bursts with a Salpeter IMF and solar metallic- 
ity. The PMAS maps of EW(Ha) em in our sample of LIRGs (see 
Figure 1), except for the nuclear region of IC 860 (Table 4), in- 
dicate ages of the young stellar populations of between 5 and 
~ 10 - 20Myr (see e.g., Leitherer et al. 1999, and Section 3.1). 
We thus chose ages of the ionizing stellar populations of between 
1 and 20Myr, and > lOOMyr for the evolved one. Both stellar 
populations were extinguished with a range of values of E(B-V). 

Combinations of evolved stellar populations with ages be- 
tween 0.7 and 5 - 10 Gyr with relatively low extinctions, and 
ionizing stellar populations with higher extinctions provided 
reasonable fits to the optical continua of most of our regions. 
The spectra of our LIRGs tend to show deeper Mglb and G- 
bands than the ULIRGs in the Rodriguez Zaurin et al. (2010a) 
sample. This suggests that the stellar populations in our sam- 
ple are somewhat older. In fact, we found little evidence for a 
strong contribution to the optical light from intermediate-aged 
stellar populations with ages of ~ 100 - 500 Myr in our sam- 
ple of LIRGs. The only exception in our sample is NGC 7771. 
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Several of the regions studied in this galaxy, including its opti- 
cal nucleus, show clear indications for an important contribution 
from an intermediate-aged stellar population (~ 100 - 700 Myr). 
Additionally the optical nucleus of NGC 7771 contains a large 
number of stars with ages in the range 30 - 50 Myr, as indicated 
by the presence of strong He i absorption features. 

The majority of the selected regions in our sample of LIRGs 
have a significant contribution to the 4600 A continuum light 
from the ionizing stellar populations, and for about one-third the 
ionizing stellar populations have a dominant contribution. The 
ages of the evolved stellar populations of our LIRGs are within 
the ranges observed in massive spiral galaxies. The fitting of the 
stellar continuum and the nebular properties yielded ages of the 
ionizing stellar populations of between 5 and 9 Myr, and extinc- 
tions in the range E(B - V) young = 0.2 - 1.8. The extinctions 
from the stellar continuum fits and the Balmer decrement agree 
with each other, and tend to be higher than those required for the 
evolved stellar populations. 

We used the brightest emission lines to study the spatially 
resolved (on a spaxel-by-spaxel basis) excitation conditions in 
the central few kiloparsecs of our sample of LIRGs. When using 
traditional diagnostic diagrams involving bright optical emission 
lines, the location of the spatially-resolved line ratios varies from 
LIRG to LIRG. Some galaxies have most of their spaxels in the 
H ii region of the diagnostic diagrams, while others have spax- 
els located both in the Hn and LINER regions of the diagrams. 
In particular, a significant fraction of spaxels with EW(Ho') em < 
20 A have enhanced [Nii]i6584/Ha, [S ii],M6717,673 I/Ha and 
[Oi]/16300/Ha line ratios, when compared to those of Hn re- 
gions. These spaxels tend to be associated with regions of dif- 
fuse Ha and Paa emission in these galaxies, that is, they do not 
coincide with high surface brightness H n regions identified in 
the HST/NICMOS Paa images. 

We also produced spatially resolved [S n]/Ll6717,6731/Ha' 
vs. [N ii]/16584/Hq' diagrams for each of the galaxies, and the full 
PMAS sample. These diagrams have the advantage of contain- 
ing more data points than the traditional diagnostic diagrams, 
and of having smaller corrections for the presence of Ha stel- 
lar absorption than diagrams using H/3. Even so, the correc- 
tions for stellar absorption to the observed line ratios are sig- 
nificant for spaxels with EW(Ha) em < 20 A. After the correc- 
tion for stellar absorption, there is still a significant number of 
spaxels in these diagrams whose line ratios cannot be explained 
as produced entirely by photoionization by young stars (e.g., 
log[S ii]/l/l6717,6731/Ha' > -0.2). These enhanced ratios are 
probably due to the combined contributions of H n region emis- 
sion and DIG emission. We estimated the fraction of DIG emis- 
sion in the central few kiloparsecs of our LIRG using two dif- 
ferent methods. The first one used Hn region photometry from 
the NICMOS Paa images. The second method was based on the 
comparison of the spatially-resolved [S n]/l6717/Ha: line ratios 
of our LIRGs, and those of the MW H n regions and DIG. The 
DIG fractions over the central few kiloparsecs vary from galaxy 
to galaxy, but are generally < 60%, as found for starburst galax- 
ies by Oey et al. (2007). 
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Fig. 1. (a) NGC 23: PMAS observed (not corrected for extinction or underlying stellar absoption) maps of emission line ratios: 
[O m]A5007/Hj3 (top middle), [O i]A6300/Ha (top right), [N n]i6584/Ha (bottom middle) and [S n]/U67 17,673 l/Ha (bottom right). 
Also shown are the maps of the observed Ha flux (top left) and EW(Ho-) em (bottom left) of the line in emission. The FoV of the 
PMAS maps is 16" x 16", and the orientation is North up, East to the left. The crosses on the maps mark the peak of the PMAS 
optical continuum emission at ~ 6200 A (see Paper I for details). The Ha maps are shown in a square root scale to emphasize the 
low surface brightness regions. 
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Fig. 1. (b) As Fig. la but for MCG +12-02-001. 
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Fig. 1. (c) As Fig. la but for UGC 1845. 
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Fig. 1. (d) As Fig. la but for NGC 2388. 

MCG+02-20-003 




arcsec arcsec arcsec 

Fig. 1. (e) As Fig. la but for MCG +02-20-003. 
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NGC5936 




Fig. 1. (g) As Fig. la but for NGC 6701. 
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Fig. 1. (h) As Fig. la but for NGC 7469. The Hfi and Ha measurements correspond to the narrow component of the lines (see 
I for details). 
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Fig. 1. (i) As Fig. la but for NGC 7591. 
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Fig.l. (j) As Fig. la but for NGC 7771, except that the FoV for this galaxy covers the central 28" x 16" (Section 2.1 for de- 
tails). The maps are shown: Ha flux (top left), EW(Ha) em (top right), [Om]A5007/H/3 (middle left), [Oi]A6300/Hff (middle right), 
[Nn]A6584/Ha (bottom left) and [S n] A/167 17, 673 1/Ha (bottom right). 
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Fig. 2. (a) Observed diagnostic diagrams for the spatially resolved measurements of NGC 23. The line ratios are not corrected for 
extinction or underlying stellar absorption. The star symbols are the nuclear region measurements given in Paper I. The line ratios 
for individual spaxels are color coded according to the value of the EW(Ho') elI1 of the line in emission: green dots EW(Har) em > 60 A, 
orange dots 20 A< EW(Ha) em < 60 A, blue dots 5 A< EW(Ha) em < 20 A, and red dots EW(Ha) em < 5 A. The solid lines are the 
"maximum starburst lines" defined by Kewley et al. (2001) from theoretical modeling, and the thick dashed lines are the empirical 
separation between AGN and H n regions of Kauffmann et al. (2003), and between Seyfert and LINERs of Kewley et al. (2006). 
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Fig. 2. (b) As Fig. 2a but for MCG + 12-02-001 . 
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Fig. 2. (c) As Fig. 2a but for UGC 1845. 
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Fig. 2. (d) As Fig. 2a but for NGC 2388. 
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Fig. 2. (e) As Fig. 2a but for MCG +02-20-003. 
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Fig. 2. (f) As Fig. 2a but for NGC 5936. 
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Fig. 2. (g) As Fig. 2a but for NGC 6701. 
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Fig. 2. (h) As Fig. 2a but for NGC 7469. The H/3 and Ha measurements correspond to the narrow component of the lines (see 
I for details). 
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Fig. 2. (i) As Fig. 2a but for NGC 7591. 
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Fig. 2. (j) As Fig. 2a but for NGC 7771. 
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Fig. 3. Spatially resolved observed [N n]/i6584/Ha vs. [S n]AA67 17,673 1/Ha diagrams for each of the galaxies in our sample. The 
individual spaxel measurements are color coded according to EW(Ha') em as in Figure 2. The spaxels corresponding to the nuclear 
region are shown as a star-like symbol. The line ratios are not corrected for underlying stellar absorption or extinction. The error bars 
are the typical uncertainties in the measured line ratios for the smallest range of EW(Ha) em , color-coded as the data points, measured 
for each galaxy. The uncertainties in the line ratios of spaxels with EW(Ha) em greater than the plotted range are always smaller. 
The solid lines are the Dopita et al. (2006) models for evolving H n regions for instantaneous star formation and 2Z metallicity in 
all cases except for the MCG +12-02-001, UGC 1845 and MCG +02-20-003 plots where we show the Z models (see Table 7 for 
the estimated abundances). The pink, dark blue and light blue lines correspond to log % parameters of -2, -4, and -6. The dashed 
lines are isochrones for ages 0.1, 0.5, 1, 2 and 3 Myr from top to bottom. The % parameter is defined as the ratio of the mass of the 
ionizing cluster to the pressure of the interstellar medium. 
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Fig. 3. Continued. For NGC 7469 and NGC 759 1 we identified the nuclear regions as the 4 central brightest spaxels at 6200 A, as 
these two galaxies contain an active nucleus and the seeing (FWHM) was larger than the spaxel size. 



